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Abstract 

An action for simplicial euclidean general relativity involving only 
left-handed fields is presented. The simplicial theory is shown to con- 
verge to continuum general relativity in the Plebanski formulation as 
the simplicial complex is refined. This contrasts with the Regge model 
for which M. Miller and Brewin have shown that the full field equa- 
tions are much more restrictive than Einstein's in the continuum limit. 
The action and field equations of the proposed model are also signif- 
icantly simpler then those of the Regge model when written directly 
in terms of their fundamental variables. 

An entirely analogous hypercubic lattice theory, which approxi- 
mates Plebanski's form of general relativity is also presented. 



1 



1 Introduction 



It has been known for some time that in general relativity (GR) the gravita- 
tional field can be represented entirely by left-handed fields, i.e. connections 
and tensors that transform only under the left-handed, or self-dual subgroup 
of the frame rotation group.[] 

The present paper presents a simplicial model of GR with an internal 
SU (2) gauge symmetry which, at least in the continuum limit, corresponds to 
the left handed frame rotation group. The gravitational field is represented by 
spin 2 and spin 1 SU(2) tensors, and SU (2) parallel propagators, associated 
with the 4-simplices, and with 2-cells and edges constructed from the 4- 
simplices, respectively. 

This is meant to provide a step in the construction of a covariant path in- 
tegral, or sum over histories, formulation of loop quantized GR. In Ashtekar's 
reformulation of classical canonical GR |[Ash86|| , [|Ash87H the canonical vari- 
ables are the left handed part of the spin connection on space and, conjugate 
to it, the densitized dreibein. The connection can thus be taken as the con- 
figuration variables, opening the door to a loop quantization of GR | |GT8(: 

m 
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In loop quantization one supposes that the 
state can be written as a power series in the spatial Wilson loops of the 
connection (which coordinatize the connections up to gauge), so the fun- 
damental excitations are loops created by the Wilson loop operators. The 
kinematics of loop quantized canonical GR requires that geometrical observ- 
ables measuring lengths |Thi96b|l , areas [|RS95| [ |AL96a|] , and volumes [ [RS95 



AL9bb 


I- 


Lew9d|, 


l'hi96a 



l'hi96a| have discrete spectra and finite, Plank scale low- 



est non-zero eigenvalues, suggesting that GR thus quantized has a natural 
UV cutoff. One would therefore expect that a path integral formulation of 
this theory would have, in addition to manifest covariance, also reasonable 
UV behaviour. 

A step toward such a path integral formulation is the construction of the 
analogous formulation in a simplicial approximation to GR. Loop quantiza- 
tion can be applied to any SU(2) spacetime lattice^ gauge theory in which 

1 In euclidean GR the frame rotation group is SO(A) which can be written as the tensor 
product SU (2) SU (2) l (in terms of fundamental representations). Left handed tensors 
transform only under the SU(2)l factor. Examples are left handed spinors and self-dual 
antisymmetric tensors, i.e. tensors a that satisfy a IJ — e 1 klo , . 

2 The lattice need not be hypercubic or regular in any way. 
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the boundary is a finite lattice and the boundary data is the connection on 
that lattice, because in this case the states can always be expressed as power 
series in Wilson loops. Moreover, for local theories of this type it has been 
shown ||Rei94|| that the evolution operator can be written as a sum over the 
worldsheets of the loop excitations,^ that path integral. All that is 

needed for a path integral formulation of loop quantized GR is a local lattice 
action for GR in terms of the left handed part of the spin connection, and 
other fields, such that the connection is the boundary data. 

Plebanski |Ple77|| found precisely such an action in the continuum. [] Here 
a simplicial lattice analogue of this action is presented. The corresponding 
path integral formulation of loop quantized simplicial GR will appear in a 
forthcoming paper. 

The present work might also be useful in numerical relativity. Unlike the 
Regge model the present simplicial model has field equations which repro- 
duce those of general relativity in the continuum limit. M. Miller [|Mil95| and 
Brewin [Bre95| have found linear combinations of the Regge equations which 
do reproduce the Einstein equations in the limit, however these do not define 
the extremum of any known action, which presents a serious obstacle to any 
simplicial path integral quantization based on the Regge model. A second 
advantage of the present model in numerical work that its field equations 
are relatively simple in terms of the fundametal variables, whereas the Regge 
equations, when written out in terms of edge lengths, are extremely com- 
plicated. Finally, a simple modification of the new simplicial lattice action 
(given at the end of the present paper) yields a hypercubic lattice action for 
GR, which might lead to a simpler and faster computer implementation than 
a simplicial model. 

In section ^| the simplicial model is presented. Field equations and bound- 
ary terms are discussed in section ^. The continuum limit is analysed in 

3 A lattice theory is "local" if the action of a region is the sum of the actions of basic 
cells (smallest subdivisions for which the action is defined) that make up the region, and 
the action of each basic cell depends only on fields living in the cell and on boundary data. 

In [Rei94] the sum over worldsheets formulation is obtained for theories whose actions 
are functions of the connection only. However, if the connection is the only boundary 
datum for each basic cell, a local action in terms of the connection only can be obtained 
by integrating out all other variables in each cell. 

4 The left handed part of the spin connection is also the boundary data of the GR action 
of Samuel [3am87|, and Jacobson and Smolin [JS87|, [JS88|. However I shall not try to 
build a lattice analogue of that action here. 
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section |j. The results of M. Miller and Brewin on the continuum limit of 
the Regge model are reviewed, and why it is that the present simplicial 
model leads to field equations approximating those of continuum GR while 
the Regge model does not is discussed. 

Finally section [5] contains some comments on the results, and also states 
the hypercubic action. Appendix |A| gives a metrical interpretation of the 
simplicial field e s , and in appendix |B| some lemmas used to establish the 
continuum limit are proved. 



2 The model 

Plebanski gave the following action for general relativity (GR) in terms of 
left-handed fields Fli77fl pjJlVMfl fl 



AF J -i^S 4 AS, (1) 

(The euclidean theory is obtained when all fields are real). This action has 
internal gauge group SU(2), with S a 2-form and an 577(2) vector (spin 1), 
F the curvature of an SU(2) connection A, and <fi a spacetime scalar and spin 
2 SU(2) tensor. The action is written in terms of the components of these 
fields in the adjoint, or 50(3), representation of SU(2). (Indices in this 
representation run over {1,2,3} and will be indicated by lowercase roman 
letters {i,j, k, I, •■■})■ is thus represented by a traceless symmetric matrix 

On non-degenerate (E fc A E fc ^ 0) solutions these fields can be expressed 
in terms of more conventional variables. X is the self-dual part of the vierbein 
wedged with it self :0 

E< = 2[e A e] + 0i = e° A e ! + ^e ijk e j A e fc , (2) 



,/3i.../3„ 



5 The definition of exterior multiplication used here is [a A 6] Q , 1 ... Qn 
a [a 1 ...a m bf3 1 ...i3 n \, where spacetime indices are labeled by lower case greek letters 
{a, /3,7, ...}. Forms are integrated according to J A a = J A ^ Ul '" Um a ul ... Um d ,n a where 
A is an m dimensional manifold, <r u are coordinates on A, the indices Ui run from 1 to 
m, and e" 1 — Um is the m dimensional Levi-Civita symbol ( e 12 - m — i and e is totally 
antisymmetric ) . 

6 Note that upstairs and downstairs 50(3) indices are the same. 
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which transforms as a spin 1 vector under SU(2)l, the left-handed subgroup 
of the frame rotation group 5*0(4) = SU(2)p ® SU(2)l, and as a scalar 
under SU(2) R . A is the self-dual (SU{2) L ) part of the spin connection, and 
4> turns out to be the left-handed Weyl curvature spinor. The non-degenerate 
solutions correspond in this way exactly to the set of solutions to Einstein's 
equations with non-degenerate spacetime metric. 

Ashtekar's canonical variables are just the purely spatial parts of A and 
£ (the dual of the spatial part of £ is the densitized triad), and, in the 
non-degenerate sector, the canonical theory derived from (|l|) is identical 
to Ashtekar's ||CD J M9T| | |Rci95|| . Since this is precisely the sector of non- 



degenerate spatial metric it is of course also equivalent to the ADM theory 
ADM62|. However, when the metric is degenerate the canonical theory dif- 



fers from Ashtekar's |[Rei95|| . Since Plebanski's theory defines an extension of 
GR to degenerate geometries, and this extension is not the only one possible, 
I will refer to this theory as Plebanski's theory. 

The simplicial model of GR presented here makes use of a somewhat in- 
tricate cellular spacetime structure known to mathematicians as the "derived 
complex" |pVlau96| . Spacetime is represented fundamentally by an orientable 
simplicial complex A. Q The derived complex is defined by subdividing each 
4-simplex of A into 10 "corner cells", each associated with a vertex of the 
simplex, as follows. A 4-simplex v has an affine structure (i.e. it is a chunk 
of a vector space) so there is a unique constant metric which makes it a unit, 
equilateral 4-simplex. Using this metric the corner cell cp of the vertex P in 
v can be defined as the closure of the set of points in v that are closer to P 
than to any other vertex of v (See Fig. |l| a)). 

cp has one vertex in the interior of z/, namely at the center C v of which 
is equidistant from all the vertices of v. f\ The other vertices of cp live on the 
boundary of v and thus in some subsimplex. Each subsimplex is equilateral, 
so the intersection of cp with a subsimplex fi is just the corner cell of P in 
fi and the vertex of cp in /i is the center of fi. It is not hard to see that cp 
is topologically a hypercube and its vertices are P, C u , and the centers of 
all the subsimplices of v incident on P. (4 1-simplices, 6 2-simplices, and 4 

7 The simplicial complex will always be assumed to be a combinatorial manifold, so it 
has every nice property that one would expect a simplicial representation of spacetime to 
have. See (gW9| for details. 

8 An equivalent definition of the center of a simplex /i is that it is the average of all 
the vertices of \i in any linear coordinates on /i. 
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a) 




b) 




P 



Figure 1: Panel a) shows the corner cell ap associated with the vertex P of a 
3-simplex. Notice that the intersection of ap with any of the triangular faces 
incident on P is itself the two dimensional corner cell of P in the triangular 
face, and that the six corners of ap are located at the centers of the simplices 
incident on P, defining, topologically, a cube. 

Panel b) shows the cell P* dual to the vertex P in a two dimensional sim- 
plicial complex. The boundaries of other dual cells are indicated by dashed 
lines. The cells and subcells of this realization of the dual complex generally 
are not flat where they meet a boundary between simplices. 

3-simplices). 

Notice that to each vertex P of a simplicial complex one can associate a 
"dual" cell P* formed by the union of all corner cells of P in the simplices 
incident on P. The complex, A*, of these dual cells is topologically dual to 
A. See Fig b). 

In the definition of the simplicial model a central role will be played by 
the 2-cells s{av) of the derived complex which are attached to the centers of 
the 4-simplices. Each is associated with a 4-simplex v and a 2-simplex a of 
v, and is a plane quadrilateral formed by the centers of u, a, and the two 

3- simplices T\ and r 2 of v that share a. s[ov) can also be thought of as the 
restriction to v of the cell a* of the dual complex dual to a. See Fig. 

s(au) will be called a "wedge" and will often be denoted by just s or ov. 

4- simplices will be denoted by u, with some additional subscripts or markings 
to distinguish different 4-simplices. 3-simplices will similarly be denoted by 
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Figure 2: Panel a) illustrates the definitions of s{av) and the edges l(vr) 
and r(r<j). In the middle lies the center of a 2-simplex a. The corners are 
the centers of all the 4-simplices u, u', ... that share a. The curve C v C Tl C v > 
connecting C v and C v > is the edge r* in A* which is dual to T\. Similarly 
\J v s(crv) is the 2-cell a* < A* dual to a. One may think of s(av) as the 
wedge of a* in v. s{av) = a* fl v. Likewise, \{yr) = r* fl v and the radial 
edge r(rcr) = cr* fl r. 

Panel b) shows the analogous structure in a 3-dimensional complex with 
a 3-simplex playing the role of z/, 2-simplices as T\ and r 2 , and a 1-simplex 
as a. 

r, and 2-simplices by <r. O-simplices, i.e. vertices, will be denoted by latin 
capitals P,Q,R, .... Finally, \i < p signifies that fi is a subcell or subsimplex 
of p which may be a simplex, a cell of the derived complex, a cell of the dual 
complex, or a complex. 

The 4-simplices will be given a uniform orientation throughout A, and 
the orientation of each wedge s(au) will be determined by the orientations 
of a and v through the requirement that a positively oriented basis on a 
concatenated with a positively oriented basis on s forms a positively oriented 
basis on v. 

In the left-handed simplicial model of GR presented here one associates 
to each wedge s{av) an SU(2) spin 1 vector e S j, which will more or less play 
the role of Plebanski's S, field.0 The role of A is played by SU(2) parallel 

9 e al , is defined to reverse sign when the orientation of a is reversed. 
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propagators along the edges of the s(av). Specifically, there is an hi e SU{2) 
for each edge l{yr) from the center of 4-simplex v to that of 3-simplex r < v, 
and there is a k r e SU(2) for each edge r(ra) from the center of 3-simplex 
r to that of 2-simplex o < r. 

Finally, a spin 2 (2) tensor tp v (represented by a symmetric, traceless 
matrix, (p^) is associated with each 4-simplex. 

The action for the model is 



ip^ plays the role of (j) lJ . 



X E esi9 s - X e si e- sj sgn(s, s)]. 



(3) 



9 l s is a measure of the curvature on s. It is a function of the SU(2) parallel 
propagators via 

ei = tr[.r gds }, (4) 

where gg s the holonomy around ds, and the Jj are 1/2 the Pauli sigma 
matrices^ gg s and ^ s may be written in terms of a rotation vector p % as 



9ds = e lp ' J = cos-^-1 + 2isin-^-p- J 
^ s = 2 sin p 



(6) 
(7) 



(p = p/\p\)- The rotation vector is essentially the curvature on the wedge 
s, and, when the holonomy is close to one, i.e. the curvature is small, 9 S 
approximates the rotation vector.^ 

sgn(s, s) = sgn(a, a) is essentially the sign of the oriented 4-volume 
spanned by the 2-simplices o and a associated with s and s: If a and a 
share only one vertex (the minimum number when both belong to the same 
4-simplex), then the orientations of a and a define an orientation for u, 
namely the orientation of the basis produced by concatenating positively ori- 
ented bases of a and a. If this orientation matches that already chosen for 
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Jl — 9 



1 

1 



1 

-1 



(5) 



11 Note that 9 S reverses sign when the orientation of s reverses because the direction of 
the boundary ds reverses, which, in turn, means that gg s — > gg s , p s — * —p s , and, finally, 
6 S — > — 6 S . 



8 



v then sgn(s, s) = 1. If it is the opposite sgn(s, s) = — 1. If a and a share 2 
or 3 vertices they lie in the same 3-plane and span no 4-volume. In this case 
sgn(s, s) = 0. 

A nice, very explicit, formula can be given for the sum in the second term 
(H). If the vertices of the 4-simplex are numbered 1, 2, 3, 4, 5, so that 12, 
13, 14, 15 form a positively oriented basis then 

Y e si esjsgn(s,s) = - Y e PQRi epsTje PQRST , (8) 

where epQR = e s ^pQR )V ), and PQR indicates the 2-simplex with positively 
ordered vertices P, Q, R. 



3 Field equations and boundary terms 

Extremization of (§) with respect to hu ViT \ is most easily carried out by 
parametrizing variations of hi via hi + Shi = hiexp[iai ■ J]. Then 



On the other hand l(u, r) C ds(a, v) when o < r, so such a variation of 
hi induces 

9dau -> 9daue iai - J . (10) 
(Here each a < r has been oriented to match dr and r to match du, with 
the effect that I is positively oriented in ds(a, u)). Thus 

50* = ^triJjj'gds] (11) 

and 

51 A = ag t ■ Y ( 12 ) 



with 



UK 



2tr[J j J l g ds ]e si = cos^-e sj + 2sin-^-[p x e s ]j (13) 



j 



2 J 2 



Extremization with respect to h^ VT ) thus requires^ 



dl 

o = MhJ—] = Y «w (*i4) 

OKI U<T 



12 If v = PQRST, t = PQiJS then X) <r<T io a ( w ) = -w PQR + w QRS - w RS p + w S pq- 
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Figure 3: The dual 2-cell a* < A* is shown with its uniform orientation 
indicated by an arrow circling in a positive sense. The routes by which the 
w aVn are parallel transported from C Un to C a to form the u aVn are indicated 
by bold lines. 

(The field equations are numbered with *s). Similarly, extremization with 
respect to fc r ( Tcr ) requires 

u<j Ul = u<j U2 (*15) 
where V\ and v-i are the two 4-simplices sharing r, and 

U UUl i = U( 1 \k r (r er )hi(y 1T )]{ i W <n/lj . (16) 

(U^-\g) is the spin 1 representation of g G SU(2)). In general u av is w av 
parallel transported from C v along the boundary ds(av) in a positive senseQ 
to C a . (See Fig. ^j. u s is thus, like w s , e multiplied by a factor which goes 
to one as the group elements hi and k r approach 1. Note that (|*15|) implies 
that all u av for a given 2-simplex a have a common value u a . 
Extremization with respect to ip l J yields 

1% = i J2 e sl e sj sgn(s, s) oc 5 {j (*17) 
(since <p l J is traceless). 



13 ^i, ^2 are numbered with index increasing in a positive sense around da*. Note that 
the definition of the orientation of s(cr, v) in terms of that of a and v defines a uniform 
orientation on a* , since the orientations of the 4-simpliccs is uniform in the complex A. 
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Finally, extremization with respect to e s implies 

or i 

1T ± = 01 '- *n¥# E e ^(s, 5) = 0. (*18) 

oe si "JU g<J , 

What about boundary terms? Suppose the simplicial complex has a 
boundary dA which doesn't cut through any 4-simplex, so it is itself a 3 
dimensional simplicial complex. No boundary term needs to be added to the 
action (El) if the connection is held fixed at that boundary, that is to say, if 
the group elements k r on the edges r(r, a) in the boundary are held fixed. 
(r(ra) in the boundary is the intersection of r < dA and the 1-cell in the 
dual, [<9A]*, of the the boundary which is dual to a). 

Another natural field to hold fixed is u a for o < dA (since it lives on the 
boundary, unlike e av and w av which live in the internal space at C u , and thus 
always off the boundary). This corresponds more closely to what is usually 
done in Regge calculus |[Reg61|| , that is, holding the lengths of the edges fixed 
on the boundary, because e au , and therefore u a , is essentially the metrical 
field variable in the present simplicial model. (See Appendix A for more on 
the relation of e av and u a to metric geometry). In this case the action must 
be modified: In the modified action the 9 S for s abutting the boundary are 
evaluated with k r replaced by 1 Vr < dA (but the definition of u a , ([16]), is 
unchanged) . 



4 The continuum limit 

The form of the action (|3p is clearly analogous to that of the Plebanski action 
([!]). Moreover, Plebanski's field equations 

' >:, A v ; x d u (19) 



6<t>v 2 
5Ip 

SIp 



D A Ei = (20) 
F i - (jfiHj = (21) 



resemble simplicial field equations ( 1*17] ), ( |* 14| ) , and ( |*18| ) respectively. We 
shall see that in the continuum limit these resemblances become exact. (The 
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simplicial field equation ( |* 1 5| ) has no continuum analog, it is an identity in 
the continuum limit I will define). 

In order to take the continuum limit of the simplicial theory I define 
below a map of continuum fields on a compact spacetime manifold M 
into simplicial fields on a simplicial decomposition A of M, which allows us to 
represent continuum field histories by simplicial ones, and a class of sequences 
{A„}^ =0 of simplicial decompositions of spacetime which become infinitely 
fine everywhere in a nice way as n — > oo. Any continuum field history 
(A, E, 0) then defines a sequence of increasingly faithful images £l n (A, E, 0) = 
(h, k, e, ip) n on the complexes A n , and corresponding evaluations I n (A, E, 0) 
of the simplicial action. 

I will confine myself to showing that a continuum limit of the simplicial 
model reproduces the Plebanski theory for continous and continously differ- 
entiable A and E. On such fields the Plebanski lagrangian, and the functional 
derivatives of the action, are continous. Unless otherwise stated A, E, and 
will be assumed to satisfy these continuity/differentiability requirements 
from now on. 

I will begin by showing that the evaluations I n of the simplicial action 
converge to the Plebanski action in the limit of infinitely fine simplicial de- 
compositions: I n (A, E,0) — > Ip(A, E,0) as n — > oo. Unless / n (A,E,0) has 
a more and more oscillatory dependence on A, E, and as n — > oo it fol- 
lows that the variations of the simplicial action due to variations 5A, #E, 
and 50 of the continuum fields converge to the corresponding variantions of 
the Plebanski action. That this is in fact the case will be verified directly. 
The Plebanski theory is thus the continuum limit of the simplicial theory in 
the weak sense that the solutions of the Plebanski theory are the continuum 
field configurations that extremize the simplicial action with respect to vari- 
ations of the continuum fields. In terms of field equations, or more precisely 
in terms of the variational derivatives of the actions - which will be called 
the "Euler-Lagrange functions" (E-L functions), one can say that A, E, and 
solve the Plebanski field equations iff the simplicial E-L functions evalu- 
ated on their image (h,k,e,ip) n = Q n {A, E,0) integrated against the image 
(5h, 5k, 5e, 5(p) n of the variations 5 A, 5E ; and 5(f) vanish as n — > oo. 

In fact a stronger result will be proved: the E-L functions vanish more 
rapidly as v — > oo on solutions to the Plebanski equations than on non- 
solutions. The simplicial E-L functions (|* 1 7|) , ( |* 14| ) , and (|*18|) on the image 
of a continuum field configuration turn out to be simply the integrals of the 
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Plebanski E-L d- forms fll9|), (pOD, and ( f2T|) over suitable <i-cells in each 4- 
simplex, modulo corrections of higher order in the radius r v of the 4-simplex 
(where the radius is determined by an arbitrarily chosen positive definite 
background metric on spacetime). As a consequence the simplicial E-L func- 
tions of the image of any continuum field configuration, whether solution 
or not, will vanish as n — > oo, simply because the cells over which the cor- 
responding Plebanski E-L forms are integrated get smaller as the complex 
is refined. Specifically, the leading term in the E-L function (|*14|) vanishes 
as rl, that in ( |*17|) as r*, and that in ( |*18| ) as r 2 v . The integrals in each 
4-simplex are sufficient to detect any non-zero components of the E-L forms 
when the simplex is sufficiently small. It therefore also follows that a con- 
tinuum field configuration is a solution of the Plebanski field equations iff all 
the E-L functions vanish more rapidly than the non-solution rate as n — > oo. 

That the simplicial field equations converge to those of the Plebanski 
theory in this stronger sense is significant. M. Miller ||Mil95|| and Brewin 



[Brc95|l have found that the analogous result does not hold for the Regge 



model |[Rcg61[ |! One can define a natural simplicial image of a metric on a 
spacetime M by covering M with a simplicial complex with geodesic edges 
and taking as the length of each edge its metric length in M. What Miller 
and Brewin found is that on the images of most solutions to Einstein's equa- 
tions (including the Kerr solution [|Mil95|| ) the Regge E-L functions do not 
vanish any more rapidly, as the complex is refined, than on non-solutions. 
On the other hand, Miller, and Brewing also showed that certain linear 
combinations of the Regge E-L functions do vanish more rapidly than the 
non-solution rate. These combinations are essentially dlR egge /dg^ v in a cell 
built out of several adjacent simplices. Their rapid vanishing guarantees that 
GR is the continuum limit of the Regge theory in the weak sense that solu- 
tions (and only solutions) to GR extremize the Regge action with respect to 
variations of the continuum metric in the limit of an infinitely fine simplicial 
complex. 

What is going on? As the simplicial complex is refined, and becomes 
much finer than the scale on which a given 5g^ v varies, the corresponding 
variations of the edge lengths, 51, approach those due to a constant Sg^, 
which are highly correlated in spacetime. Thus extremization (modulo error 

14 Brewin defines a whole class of discrete field equations which approximate GR. One 
set of field equations in this class is a sum of Regge equations. 
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terms which vanish as the complex is refined) with respect to continuum 
metric modes is a much weaker requirement than extremization with respect 
to all edge length. It leaves many modes of the edge lengths free, because 
it allows many non-zero distributions of values of the E-L functions over 
the complex. On the other hand, the fineness of the simplicial complex 
also restricts the possible modes of the edge lengths in images of continuum 
metrics. On a suffiently fine complex the metric will be nearly constant, and 
the distribution of edge lengths, and similarly values of the E-L functions, will 
be severely restricted simply by virtue of the fact that the edge lengths are the 
image of a continuum metric. The question is then whether these restrictions, 
together with the field equations that come from extremizing with respect 
to the continuum metric, are enough to imply that all the simplicial E-L 
functions vanish (faster than the non-solution rate as the complex is refined). 
In the Regge model the restrictions are not quite enough. In the simplicial 
model presented here the analogous mechanism does work. 

The result of Miller and Brewin seems a serious problem for path inte- 
gral quantizations based on the Regge model. A path integral over the edge 
lengths will not have stationary phase trajectories corresponding to the full 
set of solutions of Einsteins equations, and thus will not provide a quantiza- 
tion of general relativity. One would have to integrate over some restricted 
set of simplicial spacetimes, but to my knowledge no such restricted path 
integral has been defined. (How these considerations bear on dynamical tri- 
angulation path integrals |Wei82[ based on the Regge action is unclear to me). 
In contrast, a path integral over the fundamental variables of the simplicial 
model presented here is, at least on this count, a viable model of quantum 
general relativity. 

In the present work it is shown that the smooth solutions of the simpli- 
cial theory, i.e. solutions which are images of A, S, <fi on complexes much 
finer than the scale on which these fields vary, are in fact the solutions of 
Plebansk's field equations. What is not shown is that smooth boundary or 
initial data cannot lead to solutions which are highly crumpled, in addition 
to the smooth solutions, or that solutions that are crumpled on the simplicial 
scale approximate continuum solutions on larger scales. 

Now to the details. 

Definition 1: The map : {A, E, 0) h- > (h, k, e, ip) of continuum fields 
on M to simplicial fields on the simplicial decomposition A of M is defined 
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by 

h t = 7V-/> J (22) 

k r = vjl A - 3 (23) 

u wi = I v{C a ,x)i%(x) (24) 

J a 

<p* = 4P{C V ) (25) 

with e au defined from u a via (|TJ) and (|13"D, R 

V denotes path ordering, and v(C a ,x) = U^(Ve ljc <? ) is the parallel 
propagator of spin 1 SU (2) vectors along a straight line from x to C a (ac- 
cording to the afline structure of a). 



This definition of J1a is not the only one possible. Other maps also lead 
to equivalence of the continuum limit of the simplicial theory and the Ple- 
banski theory. For instance maps such that h, k, e, ip converge to those 
of Definition 1 as the simplicial complex is refined are viable alternatives. 
However the chosen here seems to lead to the cleanest proofs. 

Before considering the simplicial complexes to be used note that we will be 
concerned with compact spacetimes or compact pieces of spacetimes, which 
always admit finite simplicial decompositions^ 

As the sequence of progressively finer simplicial decompositions I will use 
uniformly refining sequences, which are defined as follows. 

Definition 2: {A n }^L is a uniformly refining sequence of simplicial de- 
compositions of a compact manifold M if 



1) A is a finite simplicial decomposition of M, 

2) A n+ i is a finite refinement of A n , 



and, in a fixed positive definite metric go which is constant on each v < Aq 

15 The map ([K5]) which defines w s in terms of e s is invertible except when the trace of 
the holonomy around ds vanishes. However, when the connection A is continous one may, 
by choosing a sufficiently fine simplicial complex make all holonomies around wedges close 
to 1. 

16 This follows from the arguments on p. 488 of |SW93| . 
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(in linear coordinates on v), 



3) r n , the maximum of the radii r v of the 4-simplices v < A n approaches 
zero as n — > oo, and 

4) divided by the 4- volume of i/ is uniformly bounded for all 4-simplices 
v < A n as n — > oo. 

In Appendix B it is proven that conditions 3) and 4) are independent of 
the particular metric chosen. 

The following lemma, also proven in Appendix B, will be useful 

Lemma 1: If h is a continous <i-form on a compact manifold M, {A n }^L 
is a uniformly refining sequence of simplicial decompositions of M, and go is 
a positive definite metric constant on each v < A , then Ve > N can be 
chosen sufficiently large so that for any <i-subcell c of a 4-simplex v < A^ 



where hc v is the constant ci-form (in linear coordinates on v) which agrees 
with h at C„, and r v is the go radius of v. 

Lemma 1 can be conveniently restated in terms of the characteristic tensor 
of c, which can be defined in linear coordinates x a on u by 



t c might be called the coordinate volume tensor of c, it is the <i-dimensional 
generalization of the coordinate length vector of an edge and the coordinate 
area bivector of a 2-cell. When c is a d-simplex with vertices P\, Pa+x 
t ai...a d = (p 1 p 2 )[oi...(p 1 p d4 . 1 )o«i]. Using t c equation (g§ in Lemma 1 can be 
written as 



where O(er^) denotes a quantity Q c that vanishes faster than rf, as n — > oo, 
that is, max{Q c /rf\c < v < A n } — >• as n — >• oo. 

Now we are ready to prove that the simplicial action (^|) converges to the 
Plebanski action in the continuum limit. 




(26) 




(27) 




(28) 
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Theorem 1: If {A n }^L is a uniformly refining sequence of simplicial de- 
compositions of a compact, orientable 4-manifold, M; A, E, are Plebanski 
fields on M with E and continous and A continously differential) le; and 
I n (A, E, (f) is the evaluation of the simplicial action (Q) on the simplicial fields 
(h, k, e, ip) n = Qa„(A, E, (p) defined on A n by (|22|) - (p5|), then 

lim^IniA, E, 0) = 7 P (A, E, 0). (29) 



Proof. Choose a positive definite metric go which is constant on each v < A . 
By several applications of Lemma 1 one shows that Ve > N may be chosen 
large enough so that 

IC - F{Cut p t a l u) \<erl (30) 
\e avi - Z(C u ) a(3t tf\ < erl (31) 

Vcr < z/ < A n . Thus 

In = ^ [^iafiI^s\Cu^2' t a l3t ' r s (a,u) 
v<A N a<p 

-i<^E ia/3 E i75 | a £ tftfsgn(a,a) + AI Nu }, (32) 

where the error term AI^ U is bounded by 

\AI Nu \ < ert[10\\F\\ + (10 + I| |0| |)(| |E| | + e)]\ c „ (33) 

(Here the norm | |T| | of a tensor T£± 9n is | \T\ \ = (^rX^lt^ir--^ 1 -^). 
Since F, E, and are continous on the compact manifold M they are 
bounded, so \AI Nu \ < er^n with n a finite constant. 

The ratio of rf, to the 4-volume V u of v is uniformly bounded. Let this 
bound be i? then < i?T4 implying that |A/jv!/| < skRVu- Therefore the 
sum of the errors is bounded by 

| £ AJ^I < e«W M) (34) 

^<Ajv 
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where Vm is the go volume of M, a finite number. 
Straightforward calculations show 



J2 tftl s sgn(a,a) = 30t 



(35) 




(36) 



Thus, putting everything together, we get 



E£W*i - y^a^JlS)\cjf^ s }+AI N (37) 



v 



N ■ 




I 



To prove the equivalence of the continuum limit of the simplicial model 
with Plebanski's theory at the level of field equations I show that the Euler- 
Lagrange (E-L) functions, i.e. the variational derivatives of the action with 
respect to the fields, are just the integrals of the Euler-Lagrange <i-forms of 
the Plebanski theory over certain (i-cells, modulo corrections which become 
negligible as the simplicial complex is refined. Then, if a "continuum limit 
solution" of the simplicial model is defined to be a continuum field config- 
uration, (A, S, 0) such that its simplicial images fl n (A, E, 0) = (h, k, e, ip) n 
have E-L functions that vanish faster than the volumes of the corresponding 
(i-cells as n — ► oo, the continuum limit solutions are precisely the solutions 
to the Plebanski field equations. 

The precise statement I will prove is 

Theorem 2: Under the hypothesies of Theorem 1 and the additional condi- 
tion that S is continously differentiable the following hold for any simplicial 
field history corresponding to a continuum field history via (p2|) - (|25|) : 

1) The simplicial field equation (|*15|) tr[Jik r ^-] = holds identically for 
all n. 



2) Vz/ < A„, s, 



I < v 




(39) 
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/^ + 0(ea ,40) 

+ ,4!) 

where r is the 3-simplex in v dual to I, and in the integrals the integrands 
are parallel transported to C v along straight lines. 

3) The Plebanski field equations are fully represented by the simplicial field 
equations in the sense that if the Plebanski E-L forms are not almost every- 
where zero, then on a sufficiently fine simplicial complex the simplicial E-L 
functions will not all vanish. 



itr[hiJ t 



dig 

de si 
din, 
dh, 1 



Proof: 1) follows immediately from the definition (f24|) of u au (A, £). (|3^) 
and from the corresponding Plebanski field equations ( |T9| ) and (|20| ) and 

dl n 1 

^TiJ = ^ E e wi e ffI/i s^n((T,(7) (42) 

= (E <0 ^E i7 ,)| c i £ tftfsgn(a,a) + 0(ert) (43) 

= ^AE^k^^ + O^) (44) 

MhiJi-fj?] = E^ = ^'E^/X^ ( 4 5) 

a ^ ct<t (J<r 

= (DAE i )\c va fr1? r + 0(erl). (46) 

a/3 

s((T,l/) 



Similarly (]40| ) follows from (fHf) and the identity X) 5 . <1 , t^' 13 sgn(a, a) = 30£ a/S 



via 



(9/ 1 

C-W^E^i^(^) (47) 



<9e si 30 „ „ 



= (^-^k^^+OCerjJ). (48) 

3) follows as a corrollary of 2). Condition 4) in Definition 2, which defines 
A n , prevents the 4-simplices from having zero volume, and therefore ensures 
that in each such 4-simplex v < A n the t s , t T , t u span the spaces of 2, 3, and 
4 index antisymmetric tensors at C u . I 
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Using (|9]) and the fact that under variations of A, E, <fi the variations of 
the corresponding simplicial fields h, e, <p are given by 

6^ = 5^{C V ) (49) 
5e wi = I SZt + Oierl) (50) 

J a 

8aj = J^5A l + 0{er u ) (51) 

(where the integrands are parallel transported to C u along straight lines), 
one finds that 



Corollary: The hypothesies of Theorem 1 and E continously differentiable 
imply that, under variations of A, E, and <j) 

lim 8I n (A, E, 0) = 6I P (A, E, 0). (52) 



5 Comments 

• A hypercubic lattice action for GR can be defined in complete analogy 
with the simplicial one. One defines the centers of n-cubes as the averages of 
their vertices, and the dual complex from these centers as in the simplicial 
context. In particular, the wedges s, and the edges I and r are defined by 
replacing n-simplices with n-cubes in their simplicial definitions. (See Fig. 

E3)- 

The action of the hypercubic model is 

^n = EE e si l - Ifu e si e sj s 9n(s, s)}, (53) 

v<D s<u " s,s<u 

where □ is the hypercubic lattice and v labels 4-cubes. sgn(s, s) = sgn(a, a) 
is the sign of the 4- volume spanned by the 2-cubes (squares) a and a dual to s 
and s respectively, provided a and a share one vertex. Otherwise sgn(s, s) = 
0. In particular this means that it is zero when the 2-cubes share no vertices. 

Plebanski theory is the continuum limit of the hypercubic theory ( |53|) in 
the same sense as it is the continuum limit of the simplicial theory. (The 
proofs of section ^ go through with minor adjustments). 
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Figure 4: This is the analogue of Fig. |2| for a hypercubic lattice. Panel a) 
illustrates the definitions of s(au) and the edges l(vr) and r(rer). In the 
middle lies the center of a 2-cube a. The corners are the centers of all the 
4-cubes u, i/, ... that share a. The curve C U C T1 C U > connecting C u and C v i is 
the edge t{ in the lattice □* dual to □, which is dual to T\. Similarly U u s(au) 
is the 2-cell cr* < □* dual to a. 

Panel b) shows the analogous structure in a 3-dimensional cubic lattice with 
a 3-cube playing the role of u, 2-cubes as T\ and r 2 , and a 1-cube as a. 
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• Bostrom, M. Miller and Smolin [ |BMS94|| found a hypercubic lattice 
action corresponding to the CD J action | CDJ91|| (which is closely related to 
that of Plebanski) by following the method of Regge [ |Kcg6 1| | and evaluating 
the continuum action on field configurations in which the curvature has sup- 
port only on the 2-dimensional faces of a spacetime lattice. Unfortunately, 
the continuum action is not unambigously defined on such field configura- 
tions. Nevertheless, Bostrom et. al. present an action corresponding to a 
particular disambiguation of this expression. This action, written in terms of 
the discrete curvature variable of Bostrom et. al. is formally similar to that 
one obtains from fl53| ) by eliminating e s using the field equations. However, 
their curvature is defined very differently from the curvature (8 S ) used here, 
in terms of the fundamental fields, which in their case is a discrete connection 
1-form and in the present case is the set of parallel propagators. 

• Several authors |[ReS89||, lmm94], | holl95| |, | |Zap96 | have proposed canon- 



ical lattice formulations of general relativity using lattice analogs of the 
Ashtekar variables. It would be very interesting to compare the present 
model with these canonical theories. This would require a canonical for- 
mulation of the present model, which has not yet been found. Perhaps an 
approach similar to that of | Kha95 | would work. 

• The simplicial theory presented here converges in the continuum limit 
to Plebanski's theory, which is not equivalent to Ashtekar's canonical theory 
when the spatial metric is degenerate. Thus one would expect a quantiza- 
tion of the present simplicial model to approximate a quantization of the 
constraints of [Hei95|| rather than Ashtekar's constraints. 
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A The simplicial field e s and metric geometry 



On non-degenerate solutions Plebanski's fields A, E, have metrical inter- 
pretations. Thus simplicial fields defined via (22) - (]25D should also have 
a metrical interpretation. In fact, when Sj satisfies the field equation (|19|) 
Sj A oc and the non- degeneracy requirement A E k ^ 0, it defines 
a non-degenerate cotetrad e 7 Q (unique up to SO(3)l transformations) and 
thus a non-degenerate metric^ 

9ap = e / a 5/je J /3 . (54) 

Furthermore, when A obeys the field equation D A X; = 0, the Plebanski 
action becomes the Einstein-Hilbert action of the metric ([54]) ,Q so this metric 
is the physical metric. 

On solutions of ([19]) £* = 2[e A e] + \ where 

[e A e] + IJ = \j A e J + ~/ J KL e K A e L (55) 



is the self-dual part of eAe. Therefore, by (|24j) , the simplicial variable is 
the self-dual part of the metric area bivector of a (modulo corrections that 
vanish faster than the the simplicial complex is refined): 

e aui = I ^ + O(erl) = 2a+°< + 0(er£), (56) 

J a 

with a^. J = e 1 a e J plc^P the metric area bivector, equal to the coordinate 
area bivector in metric normal coordinates. 

If normal coordinates x 1 are chosen so that a lies in a spatial hypersurface 
(x° = constant) then e aui = |e' ; jfca£ fc - exactly the normal area vector of a. 

It would be nice to give a metric interpretation of e au also away from 
the continuum limit, to make possible a direct comparison of the present 
simplicial theory with Regge calculus [|Reg61|| . This is difficult since, unlike 



in Regge calculus, the simplices of the present model are not flat. The wedges, 
s, carry curvature, 8 S , which does not generally vanish, even on solutions. 



17 A spinorial proof of this result is g iven in |CDJM9l| . A proof in the language of 50(3) 
tensors is provided i n App endix B of | Rei95| . 
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For a proof see Rci95| , Section 3. 
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However, when the holonomies gg s are all 1 the metric interpretation of 
the continuum limit extends to arbitrary simplicial complexes. The only non- 
degenerate solutions satisfying this requirement exactly are flat spacetimes, 
even though in continuum euclidean GR there are curved non-degenerate 
solutions with vanishing self-dual curvature. It seems that the discrete solu- 
tions approximating curved non-degenerate anti-self-dual solutions necessar- 
ily have some self-dual curvature. Nevertheless, it is interesting to see how 
the metric emerges even in flat solutions. 

With gg s — 1, w s = e s , so field equation (|*14|) requires 



E e — = °- ( 57 ) 



<T<T 



(57) imposes 12 independent linear constraints on the 30 e ai/ i, which imply 



just that there is a 2-form (18 components) such that 

e cri;i = ^viaptff ■ (58) 

Field equation (|* 1 7|) and the non- degeneracy condition 

e <™ ' e ^sgn(a, a) ^ (59) 



a,cr<u 



are equivalent to T, ui A T, V j oc 6^ with A ££ ^ 0. As for the continuum 
£ fields this implies £ defines a non-degenerate metric, and that e av is the 
self-dual part of the metric area bivector. 

The roles played by the field equations are worth noting. A co-tetrad that 
is constant on a simplex defines an image of the simplex in an affine 4-space 



£ with metric 5u and a fixed orthonomal basis {^0,^1,^2,^3}. ( |57| ) ensures 
that any 3-simplex can be mapped into S such that the of e UV i of its faces are 
the self-dual parts of the area bivectors. This requirement fixes the image 
(the "geometrical image") of the 3-simplex up to SU(2)r transformations. 
( |*17D and the non-degeneracy condition ensure that the 3-simplices can all 
be mapped geometrically into £ by the same co-tetrad, so they ensure that 



19 Proof, clearly the simplex can be uniquely mapped into the spatial (123) hyperplanc 
of £ such that the e av i are the spatial normal area vectors. These are equal to the self-dual 
parts of the spacetime area bivectors, which are invariant under SU (2)r transformations. 
Thus any SU(2)r transformed image of the 3-simplex fulfills the requirement. With some 
more effort one can show that these are the only allowed transformations. 
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the image 3-simplices all can fit together to form a 4-simplex. Finally, in the 



gauge in which e av = u a field equation (|*15|) ensures that the geometrical 
images of the 3-simplex faces of neighboring 4-simplices match up modulo 
an SU{2)r transformation. The only remaining field equation, ( |*18|) , simply 
requires tp v = 0. 

The field equations thus restrict the simplicial fields to ones corresponding 
to a Regge type simplicial geometry determined by edge lengths. Moreover, 
since the transformation needed to match up the geometrical images of the 
3-simplex shared by two 4-simplices is right-handed, the holonomy around a 
2-simplex of the metric compatible connection, which by definition transports 
the one image of the 3-simplex into the other, is purely right-handed. Since 
the metric compatible holonomy leaves the image of the central 2-simplex 
invariant, it can only be right-handed if it is 1. 

In a curved field configuration ( |*14 ) and ( |*15| ) contain curvature terms 
which spoil the exact metrical interpretation of e s (though it is of course 
recovered as the continuum limit is approached). 



B Lemmas for the continuum limit 

Lemma B.l: If {A n }^L is a uniformly refining sequence of simplicial de- 
compositions with repsect to one positive definite metric go which is constant 
on each 4-simplex of A , then it is a uniformly refining sequence with respect 
to any other such metric g' . 

Proof: Since all A n n > are refinements of Ao and Ao has a finite number 
of simplices it is sufficient to prove that conditions 3) and 4) of Definition 2 
hold with respect to g' Q in each 4-simplex Vq < A . Inside Vq go and g' are 
constant metrics, and since they are both positive definite they are related 
by a non-singular linear transformation. Hence there exist non-zero, finite 
constants a and b so that 

r' v < ar v (60) 
VI = bV v . (61) 

(' quantities are calculated with g' ). 

Thus 3), r n — > as n — > oo, implies r' n < ar n also approaches zero in 
this limit, that is 3), also holds with respect to g' . Furthermore, if 3c > 
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such that V v > erf, then VI = bV v > cbrf > cb/a 4 r' 4 , so 4) with respect to go 
implies 4) with respect g' . I 

Lemma A. 2: If / is a continous function on a compact manifold M and 
is a uniformly refining sequence of simplicial decompositions of 
M then Ve > N can be chosen sufficiently large that \f(x) — f(y)\ < e 
Vx, y G v < A N . 

Proof. Since the simplicial decomposition Ao is finite, it is sufficient to prove 
the lemma in each uq < Aq separately. The continuity of / implies that, 
with respect to any constant, positive definite metric go on u , there exists 
for each point x G vq an open ball B(x,r x ) of radius r x > about x such 
that \f(x)-f(y)\<e/2VyeB(x,r x ). 

Since z/ is compact it can be covered by a finite set of balls {B(x m , l/3r Xm )} 
Now let r = l/3max{r Xm }, and note that Vrr G vq B(x, r) C\vq C B(x m , r m ) H 
vq for some m G {1,...,M}. This implies that \f(x) — f(y)\ < e Vx,y G 

r) fl vq. Since N can be chosen large enough that r v < r W < n 
the requirements of the lemma are satisfied in v . I 

Corrollary (Lemma 1): If h is a continous rf-form on M and gr is a positive 
definite metric adapted to Ao, then Ve > N can be chosen sufficiently large 
so that for any rf-cell c < v < An 

\jh-jhc v \<er d v , (62) 

where hc v is the constant rf-form (according to the affine structure of v) 
which agrees with h at C v , and r u is the g Q radius of v. 

Proof. Fix on each vq normal coordinates x a of go. 

I fh - I h Cv \ < E / \Hx) ai ... ad -h(C v ) ai ... ad \\dx a K.Jx a *\. (63) 

Jc Jc ai...a d Jc 

Furthermore, since the components h aimmm0ld are continous functions, one may 
choose N large enough that \h{x) ai _ ad —h(C v ) ai _ ad \ < e(dimension M)~ d 2~ d . 
Thus 

\j h ~j h c v \ < emax{j \dx ai ...dx ad \} (64) 
< er d u . (65) 
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